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Abstract 



We study a quantum dynamical semigroup driven by a Lindblad generator with a 
deterministic Schrodinger part and a noisy Poission-timed scattering part. The dynamics 
describes the evolution of a test particle in M", n = 1,2,3, immersed in a gas, and the 
noisy scattering part is defined by the reduced effect of an individual interaction, where 
the interaction between the test particle and a single gas particle is via a repulsive point 
potential. In the limit that the mass ratio A = tends to zero and the collisions become 
more frequent as ^, we show that our dynamics ^t,x approaches a limiting dynamics 
with second order error. Working in the Heisenberg representation, for G E S(L'^(M")) 
n = 1,3 we bound the difference between ^t,\iG) and <I>^^(G) in operator norm propor- 
tional to A^. 

1 Introduction 

The study of spatial decoherence in atom optics has inspired the derivation of certain Marko- 
vian master equations as models [16] for the reduced dynamics of particle interacting with 
an environment. A Markovian approximation for the particle is made possible in part by the 
assumption that the degrees of freedom of the environment operate on a much shorter time 
scale than the particle. In this case, the individual interactions between the particle and the 
environment are effectively instantaneous with respect to the time scale of the particle. Thus 
many of the derivations of decoherence models in atom optics begin with an analysis of the 
scattering operator for the interaction between the particle and a single member of the reser- 
voir fl2[ [HI [11] . A study of models where this scattering assumption can be made more rigorous 
can be found in [3 [H [1] . 

The current article concerns a quantum Markovian dynamics simulating the evolution of a 
large test particle of mass M immersed in an inhomogeneous gas of light, high-speed particles 
of mass m. We begin with a dynamics semigroup ^t,x in the Heisenberg picture operating on 
B{L'^(R"')) and governed by an equation of the form: 



d 
di 





where $o,a(G') = G, X = P is the vector of momentum operators, p is a positive trace class 
operator on the Hilbert space of a single reservoir particle L^(M"), and Sa G i?(L^(]R")®i^^(M")) 
is the unitary scattering operator for a repulsive point interaction. The scattering operator Sa 
for two particles which interact through a point interaction is determined through center-of-mass 
coordinates by the scattering operator for a single particle in a point potential. Expressions for 
the scattering operator for a particle in a point potential can be found in [2] . Non-trivial point 
potentials only exist in dimensions 1, 2, and 3. Our parametrization of Sa by A corresponds to 
holding the strength of the interaction fixed while varying the mass ratio and a discussion can 
be found in [5]. 

Previously it has been shown in [5j that, under some norm conditions on p, there exists a 
c> such that for all A > and G G B{L^{W)) 

\\j{Tt2[{I p)Sl{G I)Sx] - G) 

- {z[V, + XV2 + ^{P, A}, G] + XiifiG) - ^{ifil), G})) II < cX'\\G\U, (1.2) 

where Vi,V2, {A)j are bounded real functions of the vector X of position operators, is a 
bounded completely positive map of the form: 

^{G) = Y, [ dkm*^^Gm^^„ (1.3) 
where m-r are functions of X, and || ■ |L„ is a weighted operator norm of the form 
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||G|U„= 11^11 + 111X1^11 + 11^1X111+ ^ (||X,P,G'|| + ||G'F,X,||)+ J2 \\\P\""G\P\ 

O^hj^d ei+e2<3 

The forms for Vi,V2,A and ip can be found in Appendix [B] and depend on p. 

In the current article, the main result is to extend fll.2p to an inequality bounding the 
difference between the semigroup ^t,x and a limiting semigroup ^fx'- 



$t,A(G') - ^UG)\U < X\t A l)tiT\\G\U, (1.4) 

o 

t,X 



for some fixed constant kt, for all A > 0, t in a time interval [0,T], and G G i?(L^(M")). 



is a semigroup satisfying an equation of the form: 

+ X{^{^1,{G))-]^{^1,{G),V{I)}). (1.5) 

Notice that we have added the second order term |A^M[A^, ■] to the generator to complete 

the square with [-P^, ■]• The dynamics ^t,\ can be constructed through a Dyson series since 
Lp is bounded and the left term of the generator written in (11.51) describes an electromagnetic 
field and is guaranteed to generate a well defined group of isometrics even for certain classes 
of unbounded Vi, V2, and A [2]. The equations (II. ip and (II. 5p formally have the Lindblad 
form [T5] except that the Hamiltonian part of the generator is unbounded in both cases due to 
the presence of P. 
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The equations (11 .ip and (11.51) should be regarded as formal, since the generators involve un- 
bounded operators (namely the momentum operators Pj). Lindblad equations with unbounded 
generators can be made rigorous using form-generator maps C.{ipi \ G; ^2) C where 1^1,1/^2 are 
in a dense domain D C L^{W) and G E B{L^{W)). A solution : 3(1^ (W)) should satisfy 

^(^i|$,(G)V'2) = /:(V'i;$*(G);V'2), 

for $o(G) = G G B{L'^{W)) and all i/^i, e D. The form £(7/^1; G; i/^a) in antilinear in and 
linear in G and 1^2, and its action is always clear from a formal expression for the unbounded 
Lindblad generator. Form generators are discussed in P, '9] and references therein. It is shown 
in Lemma \C.1\ that a rigorous interpretation of (II. ip and (II. 5p yields the existence of unique 
semigroups ^t,x, completely positive maps satisfying ^t,x{I) = xi-^) — ^- ■^o'^ 

introductory discussion of quantum dynamical semigroups see [3]. 

The dynamics $f a can be interpreted as the reduced dynamics, in the Heisenberg picture, 
for a test particle immersed in a spatially inhomogeneous gas where the effective state of the 
gas particles seen by the test particle is p. The use of the scattering matrix Sa in the noise 
term of (II. ip assumes that the gas is sparse enough so that non-trivial interactions between the 
test particle and gas particles occur individually and that interactions occur on a faster time 
scale than the kinetic motion of the test particle. If the duration of individual interactions is 
negligible compared to the time scale of the kinetic evolution of the test particle, the relevant 
information of an interaction is the translation of the "before picture" to the "after picture" 
determined by the scattering operator. 

The limiting regime A = ^ ^ 1 corresponds to a reservoir of low mass, high speed gas 
particles. By expanding our noise term in (II. ip around A = ^ = with the test particle mass 
M and the effective density matrix of a single gas particle p held fixed, the fixed momentum 
distribution p{p,p) determines growing velocities as v^p) = ^ = aa?- For the same reason, 
small A corresponds to a high temperature regime: E{p) = = ^f^g. The noise term in (II. ip 
has A dependence in the interaction frequency coefficient which grows on the order of the 
reservoir particle speeds and A dependence in the scattering operators Sa. In the limit A — > 0, 
Sa converges weakly to the identity operator. This has the interpretation of quantum tunneling, 
since a fast reservoir particle will tend to tunnel over the potential of the test particle. Without 
the frequency coefficient ^, the noise term would then tend to zero as A — 0. Intuitively, a well- 
defined limiting dynamics emerges since the increased rate of tunneling through the potential 
of the test particle is compensated by the increased frequency of opportunities for collisions 
with the high speed gas particles. 

By expanding A in the noise term around the infinite temperature regime A = only to 
first order, our limiting model neglects a description of energy relaxation. This disadvantage 
is shared by the models fi2[ |8] which are derived in similar limits involving expressions for 
scattering and small mass ratio ^ <^ 1. Since the limit A ^ involves the particles in the 
reservoir moving comparatively infinitely fast compared to the test particle, the first order 
approximation yields an environmental noise whose action is invariant of the momentum of 
the test particle. Clearly, an environmental noise term which is invariant of test particle's 
momentum can not describe the relaxation of the particle's state to the temperature of the 
reservoir. To see the invariance mathematically, let W(^q-^pj = e^i^+'^p^ be the Weyl operator 
corresponding to a space shift by q and a momentum boost by p, then the noise part of the 
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generator in (11 .Sp satisfies the momentum boost covariance: 

^(Wfo,p-)G'W-(o,p)) - ^{v^(/),l^fo„-3G'W(o,p)} = I^(*o„-)(v^(G) - G})W(o,,> (1.6) 

We conjecture that energy relaxation will be described by the dynamics generated by a second- 
order expansion in A of the noise term in fll.ip . Existence of dynamical semigroups with 
unbounded Linblad generators satisfying momentum boost covariance fll.51) is briefly discussed 
in [TO]. 

Since p is a trace class operator, it is Hilbert Schmidt, and thus the integral kernel p{xi,X2) 
in the position basis is in L^(M" x M"). Hence p{xi,X2) is small (at least in Hilbert-Schmidt 
norm) outside of some compact set. We interpret the dynamics ^t,x as valid for a finite time 
period in which, intuitively, the test particle is contained in a region where the effective state 
p is relevant. In the limiting dynamics + ^^^2 acts as a field potential and XMA acts as 

a vector potential. The presence of the A, Vi + XV2 terms and also the form of the map (f for 
the limiting dynamics differs from those yielded by a related limit (involving a scattering 
process for a massive particle in a gas) found in [51 [11] , and the difference is due to the spatial 
inhomogeneity of the gas described by p. In their setup, the limiting evolution of an observable 
(i.e. the Heisenberg representation) Gt satisfies a differential equation of the form 

^Gt = i\—,Gt\+f dkr(k)(e'''^Gte-'''^ -Gt), (1.7) 

dt ^ Jr3 ^ ' ^ ' 

— * — * 

for a specified Poisson density rate r{k) determining the rate of momentum boosts k transferred 
to the test particle. 

We can make a connection between our limiting dynamics (11. 5p and an equation of the 
form (ll.7p by considering a parametrized collection of positive trace-class operators p^ which 
will describe a spatially uniform environment in the limit e ^ 0. Define p^ through integral 
kernels in the momentum representation as 



/9,(fci,4) = e"^""'^ / dkP{k) 



e 



V2 



rce 



for some probability density P{k) on. As e — > 0, ^^(^1,^2) converges to zero for ki 7^ k2 
and asymptotically has the order of e'^"^^^^ P{k) for k = ki = k2- Using the expressions in 
Appendix [HI we see that Vi(e), V2(e), and A{e) converge to constants (and thus can be gauged 
out), and fe{G) converges to an expression as in (II. 7p which is proportional to 



[ dkPik) [ rf{;e*^(-"^+^)Ge~'^(~"^+^), (1.^ 

Jr" J\k\ = \v\ 



'\k\ = \v\ 

SO that ipo{G) — ^{lpo{I) , G} has the form of the noise in (I1.7P . The natural interpretation is that 
an incoming reservoir particle of momentum k in momentum transfers of —v + k with v on the 
shell of radius \k\. The Kraus decomposition suggested by (11.81) is an integral combination of 
conjugation by unitary Weyl operators corresponding to momentum shifts. This contrasts with 
the general situation for (p in which the terms rrij ^ determining the Kraus decomposition ( IB.Sp 



are non-unitary combinations of Weyl operators e^-^^ ^^^^ for fixed k and varying v on the 
surface \v\ = \k\. 
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The dynamics ^t,\ can be constructed as 

oo 

$,,A(G') = e-x^- / F,.,„Vx---VxF,,{G), (1.9) 

n=0 JO<ti-t„<t 

where Vx{G) = Tr2[(/ ® p)Sl{G ® I)Sx] and Ft is generated by «[^, ■]. Since Pa and Ft are 
completely positive maps, compositions and sums of them will be also, and hence the $t,A's are 
completely positive. The construction (11.91) lends to the interpretation of the dynamics ^t,\ as 
describing a free particle (no field potential) with Poisson timed collisions occurring at rate y 
with gas particles with statistical outcomes governed by V\. The limiting dynamics can 
be constructed as a Dyson series using the group F^ (defined below) and perturbative part 
X{(f{-) — ^{ip{I), •}). However, the complete positivity of the maps is not apparent without 
the analysis in the proof of Lemma IC.li 

The strategy of our proof of (11.41) is based on showing that ^t,x has second-order in A error 
for an integral equation that the dynamics solve. The basic idea is that a process close to 
solving an integral equation is close to the solution of the integral equation. We show that the 
error £'j($. a(G')) for the process ^r,x{G) at time t defined by 

i?;($,A(G)) = F^iG) + A£dsF,o_,(vp($.,A(G)) - i{$,,A(G'),^(/)}) - $,a(G), (1.10) 

where F^ is the unitary evolution generated by il^iP + XMA)'^ + Vi + AV2, ■], is uniformly 
bounded by a constant multiple of A^||G||i„„ for all G G i?(L^(M")) and times f in a finite time 
interval [0,r]. With ^t,\ replaced by in (11.101) . the error on the left side of the equation 
becomes zero, and we can apply the elementary inequality flA.ip to bound the norm of the 
difference between ^t,xiG) and ^^ ^^{G). 

A solution to a differential equation can be written as the solution to different integral 
equations and we choose (11.101) because the perturbative part y?(-) — ^{(p{I), ■} is a bounded 
map. However, using the bound (II. 2p . it is more natural to bound the error Et{^.^x{G)) in 

E,($,a(G')) = FtiG) + j\sFt^s{i[Vi + XV^ + \{P.A + yMi'2,$,,A(G')] 

+ Ay.($,,A(G)) - ^W), $.,a(G)}) - $*,a(G). (1.11) 

Since the perturbative part of the integral equation (11.111) involves unbounded terms (the 
momentum operators {P)j are unbounded), we can not apply Proposition lA.ll here. Thus we 
use Proposition IA.2I to relate the error Et{^.^\{G)) to the error E[{(^.^\{G)). A basic necessity 
of our technique is that we will need to verify that the semigroup (^t,x is minimally well behaved 
with respect to the weighted norm || ■ ||^„„ in the sense that there exists an a > such that 

\\^t,x{G)\U<e^'\\G\U. (1.12) 

for all G e S(L2(M")), t, A > 0. 

In Section [21 we develop general conditions to guarantee that a semigroup obeys an in- 
equality of form (11.121) for some weighted operator norm. Section [3] checks these conditions 
specifically for ^t,x and || ■ \\wn and gives a bound of the type (11.121) which will, importantly, be 
independent of A. Finally, Section H] contains Theorem 14. 2[ the proof of which gives the details 
for the inequality (II. 4p . 
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2 Bounding weighted norms of evolved observables 



In this section, we develop conditions for attaining inequalities of the form 

||$,(G)||Y<e"*||G||Y, (2.1) 
for a > 0, where $t : B{T-C) B{7i) is a semigroup satisfying an integral equation 

^t = Ft+ [ dsFt.s^^s, (2.2) 







-Fs,\E' : BiTi) — > B{T-C) for a complex Hilbert space Ti, Fg forms a group (with a possibly 
unbounded generator), is bounded, and || ■ ||y is an weighted operator norm determined by 
densely defined closed operators Yij, Y2J, j = 1, ■ ■ ■ , n as 

n 

\\g\\y = \\g\\ + J2\\Yi,jGy;j- 

Since the adjoints Y^j and l^j are densely defined, the operator YijGY2 j should be interpreted 
as determined by a densely defined quadratic form F{'i/ji;i/j2) = (Y^ jifjilGY-^ ^11)2) when F is 
bounded. Clearly ||G||y may be infinite for some G G B{T-C), so we look a the subspace 
A C B{T-[) with finite norm. The inequality (12. ip just shows that the dynamics $4 is minimally 
well behaved with respect to the weighted norm || ■ ||y. 

Our technique for bounding ||$t(G)||Y relies on the assumption that maps '^Y,j,s of the form 

^y,iG) = ^iY,,GY2:^) - Y,,^iG)Y2:^ (2.3) 

are bounded. In practice, it may be useful to shift this condition to certain "dilated" operators 
related to Yij and Y2J. For instance ii Ti = L^(M") and Yij = Y2J = \P\, then it becomes 
possible to look at (12.31) with Yij replaced by P and Y2J replaced by (P)*. The dilations are 

natural, since we have intertwining relations such as t^P = {P — k)Tk for = e**^"^, which we 
use in the analysis of Section[3l In this situation, would be a map from B{l-L) to BiH^C""), 
and \E' is interpreted to act component-wise. We formalize the concept of a "dilation" of an 
operator in the following definition. 

Definition 2.1. Let A be a closed operator with a dense domain in a Hilbert space Ti. A 
dilation of A is a closed operator 

A:D{A)cn^n®ndii, 

where Tidu is a Hilbert space and A satisfies 

A* A = A* A. 

The following lemma shows how dilating an operator in certain expressions leaves the norm 
unchanged. 

Lemma 2.2. Let A, B be closed operators on a Hilbert space Ti with dilations A, B correspond- 
ing to dilation spaces H^^J and Tif^i respectively. For any G G B{H), 

\\AGB*\\Bin) = \\AGB*\\^^^^^w^^^^,^^y 

where the norms can be infinite. 
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A linear map \1/ : -B(7Y) —>■ B{T-C) is said to be dilation bounded if 

\\^ <^ In II = c < oo 

1 1 ^ fTaux 1 1 

for an infinite-dimensional Hilbert space Haux- Clearly for such a map, luauJl — ^ 
when Tiaux is finite-dimensional. In particular, a bounded completely positive map is dilation 
bounded. 

Proposition 2.3. Let $f he a semigroup determined by \2.^) where ^ has a dilation hound h, 
and Ft forms a group. Assume that Ft is an isometry and that $f is an weak contraction with 
respect to the operator norm. Let || ■ ||y be a generalized operator norm and lij, l^j be dilations 
ofYij,Y2j with dilation spaces Tii^du, 'H2,dii- Suppose also that: 

1. there exists some constant a G M+ such that 

||F,(G')||v<e"*||G||Y, (2.4) 

2. and defining the commutator j{G) = Yij'^{G)Y2,j ~ i^hjC^Y^^) for j = 1, ■ ■ - n, there 
exists a c > s.t. for all j and any G 

||^y,(G)|| <c||G'||y. (2.5) 

For a = a + no, we have the following inequality: 

\\MG)\\Y<e-'\\Gy. 
Proof. Recall the Dyson series identity 

Yi,MG)Y2:, = MYiMlj) + iys<^t-s{F:^s'^y^^^y^Ft-s)<l>s{G). (2.6) 

Taking the operator norm of both sides of (12.61) . and using the fact that $t is contraction w.r.t. 
the oo-norm and that Ft is an isometry, we get the inequality 

\\Y,,MG){%,,r\\ < \\Y,,,Ft{G){%,,r\\ + l^syy^^^y^^^^^ 

< \\Y,^jFt{G)Y2*J+c f dse^'-'^^<l>,{G)\W. 

Jo 

Applying Lemma [2.21 we get 

\\Y,j^t{G)Y*J < \\Y,jFt{G)Y*J\ + c f ds e^'-''>^\\<^s{G)\\Y. 

Jo 

Now if we sum over j, add ||$f(G')|| to the left and HGH to the right (since $j is a contraction), 
and apply (12.41) . then we arrive at 

||$t(G)||Y < ||G'||Ye'^* + nc / dse^^'-'^\\<^s{G)\\Y. 

Jo 

Finally, applying Gronwall's inequality, we conclude 

||$*(G')||Y<e('^+"^)*||G||Y. 

□ 
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3 Bounding \\^t,x{G) 




The goal of this section is to apply Proposition 12.31 to ^t,\ and || ■ and to get a bound of 
||^t,A(G')||^n independent of A. ^t,\ satisfies the integral equation 

<l>t,x{G) = Ft{G) + [ t;sF,_,^A<f.,A(G'), where 

^,(0) = j{TT2[iI ® p)SUG ® I)Sx] - G) and F^iG) = e^^^'^'^G). 

In our setup for an application of (12. 3p . the weights |X| and \P\''" in our norm || ■ are 
dilated as X : L'^{W) L2(M") ® C" and P®'' : ^^(M") ^ L2(R") ® C"*" respectively, which 
are convenient for verifying condition (12.51) . To check the condition (12.51) we need to bound 
commutation maps 

• ^2,a(G') = X^a(G')-^a(XG), 

. ¥^}iG) = ispyx^xiG) - ^xiisPyxG), 

where 5" e M„(M) (and we only really need the case S = For \E'j A, j = 2,3 analogous 

maps are defined with transposed vectors of operators multiplying from the right however the 
analysis for these is the same. 

The proof of the following proposition will require notation and results from [5] . In particu- 
lar, we make use of the forms derived for Tr2[(/(S>p)A^] and Tr2[(/(8>p)A^(G®/)A] in [5], Prop. 
2.2] and [5], Prop. 2.3] respectively for Aa = Sa — /• Controlling the integrals of operators in 
these derived expressions requires the use of [5, Prop. 3.1] and [5", Prop. 3.3]. 

Define the weighted trace norm 

\\p\Un = Ml + J2 E iii^r-'+i^.,[^„p]]iii 

n 

+ 5^5^|||P|"-2+e^^.pX,|Pp-2+e||^^|||p|2(„-2)^|P|2(n-2)||^^ (3.1) 
e j = l 

where the sums in e are over {0,1} and { — 1,0,1} in the one- and three-dimensional cases 
respectively. 

Proposition 3.1. Let p satisfy that \\p\\wtn < oo, then there exists a > such that for all 
0<X,Ge fi(L2(R")), 

||<fu(G')IUn<e"lG|U„. 

The norm || ■ \\wtn is somewhat stronger than what is required for the proof, but it does not 
help us to improve it since ||p||u,tn is used in our application of the result from [5J. 

Proof. By the construction (II. 9p . the maps are ^t,x completely positive. Also from (II. 9p we see 
that ^t,\{I) = I- Hence the maps $j_a are contractive since ||$t,A|| = '^'t,A(-^)|| = ||-^|| = 1- 

Since Ft leaves P invariant and Ft{Xj) = Xj + tPj, using the triangle inequality and that 
1 + t < e*, it follows that there exists an a > s.t. ||Ft((j')||^„ < e'^*||G||i„„. Thus condition (1) 
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of Proposition 12.31 holds. Checking the second condition (2) of Proposition 12.31 requires many 
calculations where the basic strategy and techniques are similar in the different cases, so we 
focus on ^1 for s = 0. 

^A,i(G') can be written as ^i,x{G) = ^l.aIC) + ^c,xiG), where 

^l,a(G') = (P«''Tr2[(/®p)A*] -Tr2[(/®p)A*]P®'')G', and 

^^cAG) = P«''Tr2[(/ ® p)AliG ® I)Ax] - Tr2[(/ ® p)Al(P«''G' ® I)Ax]. 

The above terms "^lA^) "^cxiG) are naturally handled individually. Using [5, Prop. 2.2] 
and [HI Prop. 2.3] the ^l,a(G') and \1/c,a(G) terms can be written as 



^ J MP JSOn 

S{\k\)[{P-{a- I)kf^ - P^'^jG 

and 



^cAG) = Y.\ [ [ da, da, Ut ^ m^. 



'SOnXSOn 

respectively, where Wj:^)^{P) = ^i^P — j^k. 

Starting with the ^E'l,a(G') term, (P — {a — I)k)'^'~ — P®'' can be represented as a sum of 
operators X^^ ® ■ ■ ■ ® X^^, where Xi = k — ak and X, = P, and at a least a single ej is 1. 

II^l,aG1| < E II / dk [ daT^,Ti^^p{d^,AP)^^K.AP) + - m) 

J'^" J SOn 



ei,--e 



jS{\k\)X,,0...X,,^G\\ (3.2) 

For each X^^ ■ ■ ■ X^^ we can unitarily rearrange all X, = P terms to the front of the tensor 
product since this is merely equivalent to a formal rearrangement of the tensors in the product 
L^(R")(8)(C")'^'^. Focusing on a single product with p < r components of P, it can be rearranged: 

dk / darrr*. ' ' \ p(ar ,(P),d^. ,(P) + (a - I)k) 

i?i(P,^,l,A)P«'(/+|P|)G'||, (3.3) 

where we have rewritten —ak + k = d^^{P) — {d^^{P) + ak — k) to match the matrix entries 
of p. Now to set things up for an application of [5, Prop. 3.1]: 

ri = (V3+|P|)"-'|[P,mp)|, 
nr = Ei(P,k,l,X) K''^' "^^^ , 

J[P,-Y-n_^P)\ 

%a = '^k,aM^k,xiP)^Hax{P) + {a-I)k), 



where n^,^^ is treated as a map n^^^^_^ : L2(M") oC"^ L2(M") (g)C"'' acting on P'^" {1 + \P\)G. 
Hence (13 21) is bounded by a constant multiple of 

||r/||i|||P|P(/+ |P|)G||. 

For ^c,x{G), we need to pick a way of rewriting the expression {'w^^_^ ^(P))®"^ — P®'^. Using 
the relation Wj:^^{P) = P — X'^^^^P — ^^^^^^^^^^ o- a(-^) ' '^^'^ write 

where Xi = P, X2 = — A^^^^P, X3 = ■^^-'-^qf^— ^"^g.^ a(-^)' ^'^'^ ^'^^ ^^'^ cases where 

not all operators in the product are P. 

ycAG)\\< J2 dk[ da,da,Ul^^^m^^^,^^^AP)* 

/ , N ^ J SOnXSOn ' ' 



We break the analysis into two cases: = 3 for some i and 7^ 3 for all i. 

For the case in which = 3 for some i, it follows that the net power of P in the product X^^ ® 
• • • ®X^^ is strictly less than r. Formally rearranging the order of the tensor product L^(]R") 
(C")®^ so that the Xi components come first, X2 components second, and X3 components last: 
r-mi -m2 > 0. Define the operator P„,,,,„a G B{L^iW) ® (C")®'"^^'"^), 
which operates on a simple tensor product f{z) (g) (g) ■ ■ ■ as 

((J -^)^ 

B.ra2,m^,xf{z) ®Vi®---Vm = f{z) ® — ^1 ® — ^ma ® ^ma+l ® " " " ^^mi+ma • 

a + X a + X 

We can rearrange our expression as: 

E^iP, k, ai, 1, A) P™,,^,,, P-^+"^ (1 + |P) G 5A(|rf^,,(P) I) "^,,,V„a(^) ^feV„All ■ 
Hence we can apply [5", Prop. 3.3] as 

v^-\k) = {Sn,3 + \kr')\kr"'^-'^'f,{k), 



:^.r — m-\ —mo 
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where the integral is acting on p®'"^(^/ -|- IPQC Hence this expression is bounded by some 
constant times 

||(5n,3 + |Pr-=^)|P|'-- 1— V|^r""^"™H'5n,3 + | ^1""') || 1 1| | PP^^ + 1^1)^11. 

This concludes the case when at least one of the is 3. Now we can move on to the 
somewhat simpler case when ej 7^ 3 for all j. Since at least one of the ej is 1, we have a 
non-zero power A™^ coming from the X2 = — A^^^4^P terms. The factor j in front of the 
expression is then cancelled, and it has the form 



5^A--^ [ dk [ da,da,Ul^^^m^^^,^^^,{P)*U\du^^^^^^ 

In this case we can apply [3 Prop. 3.3] as 

nfljP) = A™-im^.,:,^^,(P)5,(M";:,(P)|)S„^^^^,,, 



where the integral of operators is acting on P^^'G, and we get that the above expression is 
bounded by some constant multiplied by ||p||i||P®'^G'|| = |||P|''G'||. 

□ 



4 The Main Theorem 

Before we get to proving the main result we need to state two inequalities similar to (11.21) . 
except that the errors is of first order rather than second order. Their proofs are simplified 
versions of the proof of (11.21) so we omit them. 
Define the weighted operator norms 

• l|G||a=||G|| + |||P|G|| + ||G|P|||, 

• ||G||b=||G|| + |||PpG'|| + ||G|P|||, 
and the weighted trace norms 

• iipIIm = iipiii + \\\x\p\x\\\, + E,E,- iii^r-'+i^.,p]iii + iiipr-vi^r-'iii, 

• iipiib,i = iipiii + EeE,, \\\pr'[p^, [x,,p]]\\i + iiip|"-vipr-i)iii, 
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where the sums in e are over {0} for dimension one and { — 1,0} for dimension three. Note 
that if we represent the integral kernel of p in the momentum basis, then IIpT^pIpT^Hi = 
^■^^dk\k\'^^ p(k^k) and [X^, ^2) = ^p{^i + h,k2 + /i)U=o- Hence the norms ||p||a and 
||p||b are like weighted Sobolev 1-norms on p. 
Define the maps "^'^ and on 3(1"^ (W)) 

Lemma 4.1. Let ^t,x be defined by 

There exists a c > s.t. for all p, G,and < A such that 

II^a(G) - < cA||p|U,i||G|U and \\^'^{G) - ^TiG)\\ < cA||p||b,i||G||b. 

Theorem 4.2 (Main Theorem). Let p satisfy that \\p\\wtn, ||p||b,i < 00. Then over any interval 
[0,T] there exists a constant kt s.t. for all t e [0,r], < A, and G e B{L'^{W)), 

\\^t,x{G) - K,iG)\\ < (t A 1)ktX^G\U. 

Proof. The basic strategy of this proof was discussed in the introduction, and now we fill in 
the details. By our conditions on p, it follows from Proposition 13.11 that there exists an a > 
such that: 

||$t,A(G)|U„<e'^iG|U„ (4.1) 

for all < A, G G B{L\W)). 

^t,x{G) satisfies the last hit integral equation: 

<!>tAG) = FtG + [ dsFt.s'^x'^sAG). 
Jo 

By adding and subtracting ds Ft^si'^x + iMX^[A^, ■])(^s,xiG)) from the above equation and 
rearranging terms, we can write: 

Ft{G)+ [\sFt,,{^l + iMX'[A',-]){<l>sAG)) - <^tAG) = Et{<l>r,x{G)), where (4.2) 
Jo 



Et{<^r,x{G)) = - f dsFt^,{^x-n)i'^sAG)) + X^ f dsFt_,i[A\.^,AG)]. (4.3) 
Jo Jo 

Eti^^.^G)) is understood as the error for the process (^r,xiG) to the equation (14. 2p . The 
process $^;^(G') solves the left side of 04.21) exactly (i.e. with a zero error). As discussed in 
the introduction, we would have liked to bound the difference between (^r,x{,G) and $^;^(G') 
by bounding the error Eti^^.^G)) so that we can apply Proposition lA.ll However, since ^\ 
contains an unbounded part involving commutations with momentum operators, we cannot 
apply Proposition lA.ll so directly. The main technical details of this proof concern using 
Proposition IA.2I to relate the error Et{(^srx{G)) to the error £'j($._a(G')) determined by the 
integral equation: 



E[{^.AG)) = F:{G)+ f dsFi^L{^, 

Jo 
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where L{G) = (p{G)-^{G, ip{I)} and is the group generated by i[jj{P+XMAy + Vi + \V2, ■]. 
The commutators with Vi and V2 could have been left with the perturbative part of the integral 
equation, but, by convention, we group them with the evolution part of the equation. 
Now we prepare a launch of Proposition IA.2[ Define the maps 

1 MX"^ 
^A,2(G) = i[Vi + X-{P,A} + XV2 + ^^,G], and 

KAG) = X{^{G)-^MI),G}). 

By Proposition IA.2t if '^'^^^Fs_r'^'x j^\^r{G), hj = 2,3 make sense and are uniformly 
bounded for s, r G [0, t], then we have the following integral relation between the errors as: 

i?^(<l>.,,(G)) = m.AG)) + fdsFl^^i^l,E,{<^.,,{G)). (4.4) 

Jo 

However, \\'^\^Ft_s'^'^^j^s,x{G)\\ < c||$s,a(G') lUn for some c > since the only unbounded 
part in the expressions are the powers of Pi, which get up to second-order in the products 
'^iFt-s'^j with i,j > 2, and || • includes such terms. 

Now we begin the analysis of bounding £'j(<I>.^a(G')). The £'(($. ^a(G')) term can be bounded 

as: 

wm-Aom < fdsm,-n){^,AGm + >^'M fds\\[A',^sAG)]\\ 

Jo Jo 

< X' f ds{Ci\\<^sAG)\U + C2\\<^sAG)\U) < A^(t A1)C;||G'|U„, (4.5) 
Jo 

where the second inequality follows from (11. 2p and the fact that A"^ is a bounded operator for 
some constants Ci, C2. The third inequality comes from (14. ip and of course that HGH < HGII^n. 
It follows that £'t($. a(G)) is uniformly bounded and second-order in A over finite time intervals 
for some constant C[ and is linearly small in t for times near zero. 
The other expression on the right side of (14. 5 p has norm less than 

fds\\n,2Es{<^-AG))l 
Jo 

It is not immediately clear how to bound this quantity since 2 involves unbounded terms 
involving the vector of momentum operators. We will reorganize these terms to a form that we 
can analyze. Using c.c.r. {P, A} = 2AP+^j A- or {P, A} = 2PA+^. A- where A- = [Pj, Aj]. 
Using these equalities, we can rewrite \I'^ 2 the following: 

1-1 1 " 

nAG) = m + A^2 + -X'A' + A- ^ 4, G] + -XJ2 KA^PjG - GP,A,), 

j i=i 

where Vi + XV2 + ^X^A"^ + A| Aj is uniformly bounded in operator norm for bounded A > 0. 
Define 

n 

b = 2m\\ + 11^211 + 111^11 + ^ II4II) A sup \\A^\\. 
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Then 



*rf.||vl/l2(i?*($.,A(G')))l| <fe rrf.l|i?.($,A(G'))l| 







n 
.1=1 



*rf.(||P,i?,($,,(G))|| + ||i?,($,,(G'))P,||). 



We have aheady analyzed Et{^.^\{G)) and the first term is bounded by C^'tHGH^n. Now let us 
look at the expression \\PjEt{^\^r{G))\\. It can be bounded by 

< \\P E,{<^.AGm 

< ds IIP (VP, _ vI;-,)($,,,(G)) II + y M £ ds IIP, [\A\\ M/,,,(G')] || . 

The second expression on the right can be bounded in operator norm by a constant multiple 
of A^tjlGjl^n, since the 1^4^ operator is bounded and the norm || ■ includes weights by the 
momentum operators Pj. The first expression on the right side takes a little more work. Writing 

P^x{<^t,x{G)) = ^x{P<^xAG)) + ^'A($t,A(G')), and 

Pni'^.AG)) = mP^tAG)) + ^Ti^tAG)), 
by Lemma [4.11 there exists a C3 s.t. 

- ^l)P<!>rAG)\\ < C3A||P$,.,(G')|U < \Gs\\<^rAG))\U < \Gse"'\\G\U, (4.6) 

and there exists a C4 s.t. 

||($; - ^';^)$,,,(G)|| < C4A||<l>,,A(G')||b < C4A||<1>,,a(G')|U„ < AC4e°n|G|Un. (4.7) 

So we have 

||P,E,(<1>.,,(G))|| < r drXiC, + C,)\\^r,xiG)\U- 



By similar reasoning, the terms ||i?s($.,A(G'))-ft'j|| have the same bound. So from Proposi- 
tion (14.41) . we have found that there is a constant C[' such that 

sup iii?:(<i>.,A(G))ii < A2(tAi)criiGiu.. 

s€[0,t] 

Hence ^t,\{G) has a small error term for the integral equation with free evolution P* and 
perturbation ^E'a 3 • 3 bounded so we can apply flA.ip to get a bound for the distance 
between $j a(G') and the solution $^,(G) of the integral equation. 

sup ||$,,a(G)-<1>:,,(G)|| <e*ll*^>3ll sup ||P2(<f ,a(G'))|| < t/t^A^HGIU^, 

se[0,i] s£[0,t] 

where Kt = e^li'^^-^UCf. Thus for any finite interval [0,T], \\^tAG)-^lAG)\\ < (tAl)/tTA2||G'|U„. 

□ 

APPENDIX 
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A Inequalities involving errors to integral equations 

Th basic idea of the following proposition is that a process having small error in solving an 
integral equation is close in norm to the solution of the integral equation. 

Proposition A.l. Let ip, Ft he hounded linear maps on BiTi) and Ft he he a semigroup with 
II FJ = 1. Let H{s) he a process H : M"*" — >• B{7i) that is bounded in norm uniformly on finite 
intervals. Define the error 

Et{H) = FtHiO) + f dsFt_MH{s)) - H{t), 
Jo 

and let H{t) be the solution of the above with Et{H) — 0. Then 

sup \\H{s)-H{s)\\ <e*ll^ll sup 11^.(^^)11. 

se[o,t] se[o,t] 

Proof. By the iterative definition H^'+^it) = Ft(i7"(0)) + ds Ft-sH''{s) with H\t) = H{t), 
jjn+i _^ jj i^niformly over finite intervals. By a telescoping series occurring on the left side of 
the following equation we have: 

Ft{Eo{H)) + J2 I dti---dtnFt-tr.^---ifFt,{Et,{H)) 

J0<tl< - <tn<t 

= Ft{H{0)) + Y, [ dtr---dtn Ft.t^ip ■ ■ ■ ipFt, {H{0)) - Hit) = H{t) - H{t). 

^^-^ Jo<ti< -<t„<t 

Taking the norm of both sides of the above equation, then 

e« sup \\E,{H)\\>\\H{t)-H{t)\\, 

0<s<t 

which implies the result. 

□ 

Solutions to two different integral equations can represent solutions to the same differential 
equation. However, if a process H{t) is not a solution to these integral equations, then it will 
have different errors for the different integral equations. The following proposition gives an 
identity relating these errors. 

Proposition A. 2 (Integral equation error identities). Let \E'i,^'2, and \E'3 he linear maps acting 
on dense subspaces of B{7i). Also let Fi{t), F2{t) be semigroups with \\Fi{t)\\ < 1 generated by 
^1 and ^1 + ^^2; respectively. Finally, let {G{s)) e be a process satisfying: 

1. G{t) e ni>2 D{-^iFs) n n^j>2 D{-^iFs-^j) for all t and s, and 

2. G{t) j^^i^Git)) , and '^iFt^s'^j{G{s)) for i,j>2 are uniformly bounded in operator norm 
for t,s in finite time intervals. 
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Define errors E}{G), E^{G) for the process G{s) at time s as: 

Ei'\G) = Fi'^Go+ fdsFi^l{^2 + ^3)G{s)-G{t), 

Jo 



eI'\G) = Ff)Go+ / dsFl'l^sGis)-Git). 

Jo 

Then we have the integral identity: 

E^iG) = Ei'\G) + fdsF}ll<if,Ei'\G), (A.l) 

and inversely 

e['\G) = E?\G) - f dsFil^,Ef\G). (A.2) 

Jo 

Proof. ( lA.ll) follows from a little algebra involving two applications of the identity: 

Fi^\G) = F«(G) + tdsFBl^2F^'\G), (A.3) 
Jo 

and (IA.2I) is similar. 



□ 



B Expressions for Vi, V2, A, and (f 

The forms for Vi, V2 , A, and ip can be written in terms of an integral kernel for p as the 
following: 



Vi = Cn / dk\k\-^ [ dvidv2p{vi,V2)e'^^'^'-'^'\ (B.l) 

Jr+ J \iJi\ = \i^2\=k 

V2 = f dk\k\-' [ dVidV2{vi + V2)VTp{vi,V2)e'^^'^^-'^'\ (B.2) 

A = c^ f dk\k\-^ [ dvidv2VTp{k,v)e'^^'^-^\ (B.3) 

^{G) = el f dk\k\-^ [ dvidv2p{vi,V2)e'^^-'^^+^^Ge-'^'^-'^^+^\ (B.4) 

where the constant Cn depends on the dimension rz, = 1,3 and can be can found in [5]. The 
expressions are correct for dimension n = 2 also except that there is an extra term missing 
in the expression for V2- In dimension n = 1, the integrals over surfaces \vi\ = \v2\ = k are 
replaced by sums. A Kraus decomposition for 99 can be found, for instance, by taking the 
singular value decomposition p = ^j\fj){fj\ where the /^'s are orthonormal. The 
operator of (11.31) can then be taken as 



3,k 



= V^cn / dvf,{v) e^^(-^+>^. (B.5) 

J\v\ = \k\ 
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C Existence and uniqueness of dynamical semigroups 



Treating the equations (11.11) and (11. 5p rigorously is not difficult, since the unbounded portion 
of the generators are merely the kinetic part involving P. In particular, the assumption that 
||p|Utn < oo, which we have used to bound the error between ^^^x and imphes that Vi, V2, 
A, and are bounded. 

Lemma C.l. Let D = {ip e L2(M") 1 1| |P|2^/,||2 < 00} and define 

+ {HG{-^T^P' + y/)^2> + (^i|Tr2[(/ ® p)S* (G ® /)S,]V^2>, anc? (C.l) 



+ {^i\G{ - i^{P + XMAf - iV, - i\V2 - \v{I))^2) + A(^i|v.(G)^2>. (C.2) 
There exist unique, conservative, semigroups of maps (^t,\ and ■^'^^c/i that 

|(^i|$,,,(G)^2) = £(^i;$i,A(G);V'2), and |(V;i|$^_,(G)^2) = £^(^1; $^,,(G); ^2). 

Proof. The arguments for x and ^^x similar, so we will just prove the result for ^^x^ 
since it has slightly more detail. First we argue that in the domain D, the operator 

is maximal accretive. By [HI Lem. B.l] the operators Vi, V2, A, and [Pj,j4j] for j = 

1, ■ ■ ■ ,n are bounded. The only unbounded terms are B = i^P"^ + {A -P}; which is contained 
only in the imaginary part of L, so we just need show that B is self-adjoint with domain 
D. However, {A, P} is a relatively bounded perturbation of 21? -P^- -^y rewriting PA = 
T.j[Pj^ A] + ^P^ we have that 

ii{AP}^iii< ii^ip^^A.-iiiiiv^ib+iiiinniiFi^ii^ 

j 

for tfj E D. Applying basic inequalities |||F|'0||2 < "^ll'^lls + "^ll-P^'0l|2- So the coefficient 

in front of ||2^P^'0||2 in a relative bound inequality can be made < 1 and by [T3], S is a 
self-adjoint operator. 

The operators djAj, Vi, V2, A, and </?(/) are all bounded when ||p||^in < 00. By [9j, it 
follows that a process Gt G -B(L^(M")) solving 

f^{iJi\GtiJ2) = C%ij,-Gucl>2) (C.3) 
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is equivalent to solving 

However, since e~*^ and Lp are bounded maps, the solution Gt to the later can be constructed as 
through a Dyson series in Wt{-) = e~*^*(-)e~*^ and operating on Gq = G. Call Gt = T^)^{G), 
then the complete positivity of the maps T^^ follows from the complete positivity of Wt and (p. 

The last thing to show is that ^^^^ is conservative. However, a Dyson series expansion 
in and ip{-) — ^{ip{I),-} solving the integral equation (11.101) (with zero error term) will 
also be a solution to (1C.3P and hence be equal to ^^^^{G). Finally, since F^{I) = I and 
ip{I) - IWil), 1} = 0, it follows that $^^(/) = / . 

□ 
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